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Abstract 

Comparing quantities to analyze charged fluctuations in heavy ion experiments the dis- 
persion of the charges in a central rapidity box was found to be best suited. Various energies 
and different nuclear sizes are considered in an explicit Dual-Parton-Model calculation us- 
ing the DPMJET code and a randomized modification to simulated charge equilibrium. For 
large enough detection regions charged particle fluctuations can provide a signal of the basic 
dynamics of heavy ion processes. 

1 Charge Fluctuations in Fixed Target Hadron-Hadron Experiments 

Let us look for a moment back to the analysis of purely hadronic multi-particle production. At 
fixed target experiments it was possible to measure the charges of all forward particles. In this 
way significant results could be obtained with low energies available at the seventies H, ||, ^, ^ |^ : 

• The charge fluctuations found involve mostly a restricted rapidity range. 

• Good agreement was obtained with cluster models. 

The Quigg-Thomas relation |^ Q (assuming neutral clusters) for fluctuation across a rapidity y boundary 

< 5Qly > = < (0>,- < 0>, >f >= c . diV— (1) 

was found to be roughly satisfled §, |, 0, |ll], |l|, |l|, 0, |l|, |6|, 0. To illustrate, how early this 
relation was usefull we consider 24 GeV proton proton scattering data. Taking the rapidity of the 
forward-backward border as variable the data could be presented as in flgure 1 |^. The top and bottom 
lines and points correspond to with and without correcting for the leading charges. The lines correspond 
to suitably normalised spectrum. 

The agreement could be improved when string type qq -charge exchanges between the clusters were 
added |^]. Such exchanges appear in a large class of models. Using the Dual Parton Model code 
DPMJET we re-checked the old results. For pp-scattering at laboratory energies of 205 GeV good 
agreement was still obtained. The Quigg-Thomas relation is satisfled with c = 0.70 comparing to the 
experimentally preferred value: c = 0.72. 

2 Charge Fluctuations in Heavy-Ion Scattering Experiments 

In heavy ion scattering it is a central question whether the charges are distributed just randomly or 
whether there is some of the initial dynamics left influencing the global flow of quantum numbers. 
The charge flow measurements could again be decisive. It is not an impractical conjecture. In 
heavy ion experiments the charge distribution of the particle contained in a central box with a 
given rapidity range [— ?/max., +?/max.] as shown in figure 2 can be measured and the dispersion of 
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Figure 1: The measured dispersion (+) and the produced charge dispersion (*), which is corrected 
for leading charge flow, is compared with the suitably normalized negative (produced) particle 
spectrum at 24 GeV/c. 
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Figure 2: Kinematic region of the "central box" 



this distribution < bQ"^ > can be obtained to sufficient accuracy. For sufficiently large gaps this 
quantity contains information about long range charge flow. In comparison to the fluctuations in 
the forward backward charge distributions the charge distribution into a central box (having two 
borders) can be expected to require roughly twice the rapidity range. 

3 Equilibrium Expectations 

Within the framework of equilibrium models it was proposed to use the quantity to distinguish 
between particles emerging from an equilibrium quark-gluon gas or from an equilibrium hadron 
gas Ipls] , |l9| , po| . For a small enough box in a central region at high energies where average charge 
flow can be ignored, the (essentially) poisson distributed hadron gas yields a simple relation 

< SQ^ > = < TVchargcd > . (2) 

for any thermalised particles with charges and ±1. The inclusion of resonances reduces hadron 
gas prediction by a signiflcant factor taken |l9| , ^ to be around 0.7. 

It is argued in the cited papers that this relation would change in a quark gluon gas to 

< >^ q2 ^ >^ Q^Q ^ iVchargod > (3) 

i 

where qi are the charges of the various quark species and where again a central region is considered. 
The coefficient on the right was calculated with suitable assumptions. A largely empirical final 
charged multiplicity iVchargod = |(-/Vgiuc + 1.27Vquark + l-27Vantiquark) was used. 

It should be pointed out that the estimate is not without theoretical problems | p2| , p3[ . There 
is also a number of systematic uncertainties in the above comparison. As explained below in a 
simple approximation the result strongly depends on what one takes as primordial particles and 
how the extra quarks needed for hadronization are modelled. Considering these uncertainties we 
follow the conclusion of Fialkowski's papers |Q that a clear cut distinction between the hadron- 
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and the quark gluon gas is rather unlikely. This does not eliminate the interest in the dispersion 
as a measure of equilibration. 



4 Various Measures for Charge Fluctuations 



For the analysis of the charge structure several quantities were discussed in the recent literature. 
Besides the classic charge dispersion 

< SQ^ >=< (Q- < Q >)' > (4) 

it was proposed to just measure the mean standard deviation of the ratio R of positive to negative 
particles 

N+ TV, ^ ^ 



or the quantity F 

where Q = — is the charge in the box. The motivation for choosing these ratios was to 
reduce the dependence of multiplicity fluctuations caused by the event structure. 



5 Evaluation of the measures 



In the region of interest for large nuclei at high energies and strong centrality the charge compo- 
nent of the fluctuations dominates. In these region all measures are simply connected by the following 
relations ||l9|l : 

< A^chargcd >< 5R^ >= 4 < TVchargcd >< SF"^ >= 4 • ■ (7) 

^ ^charged ^ 

and the question of the optimal quantity is somewhat esoteric. To show this statement all three 
quantities were calculated in the Dual Parton model implementation DPMJET j25| shown in 
flgure 3. 

For the most central 5% Pb-Pb scattering at LHC energies (^/s = 6000 A GeV) there is indeed a 
perfect agreement between all three quantities as shown in the figure. This agreement stays true 
for analogous Pb-Pb data at RHIC energies (y^ = 200 A GeV) . 

Outside the region of interest - i.e. in the region of lower particle densities - the conventional 
dispersion, < SQ >, has clear advantages. The alternatives are not suitable for small Ay boxes in 
less dense events, 

• as if no particle in the corresponding box exists in rare events 0/0 or cx) is undefined and 

• as if one somehow fixes the problem (e.g. by not considering problematic events) their mutual 
relation is destroyed. For the minimum bias S-S scattering at these energies the agreement 
is lost and the new measures behave rather erratic . The same erratic behavior for the 
new measures is found for the proton-proton case (figure 4). 

As any conclusion will have to depend on a comparison of central processes with minimum bias 
and proton-proton events, there is a clear advantage to stick to the dispersion of the net charge 
distribution < SQ"^ >[] . 

^ New RHIC data of the PHENIX collaboration which appeared after the talk resp. paper|p^ confirm the 
problem with < SR^ > which does not appear for < 5Q^ >. Equation 7 does not hold in their case as a restricted 
azimuthal range was considered effectively reducing the density. For the measured very narrow rapidity range the 
data are essentially consistent with (hadronic) statistical fluctuations. 
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Figure 3: Charge fluctuations for the most central 5% Pb-Pb scattering at RHIC energies (^/s = 
200 A GeV and at LHC energies (^/s = 6000 A GeV). Also shown are corresponding data for p-p 
scattering 




Figure 4: Charge fluctuations for minimum bias pp scattering at SPS, RHIC and LHC energies 
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Figure 5: Example of a quark line graph 

6 A Simple Relation between the Quark Line Structure and Fluctua- 
tions in the Charge Flow 

To visualize the meaning of charge flow measurements it is helpful to introduce a general factor- 
ization hypothesis. It postulates that the light flavor structure of an arbitrary hadronic amplitude 
can be described simply by an overall factor, in which the contribution from individual quark lines 
factorize. It is for most purposes (which consider long range fluctuations) an adequate approxi- 
mation. 

The hypothesis can be used to obtain the following generalization of the Quigg- Thomas rela- 
tion ||2^, |8|. It states that the correlation of the charges Q{yi) and Q{y 2)1 which are exchanged 
during the scattering process across two kinematic boundaries, is just 

< SQ{yi) ■ SQ{y2) >= ncommon lines < Sq^ > . (8) 

• where the charges SQ{yi) = Q{yi)— < Qivi) > were exchanged across two kinematic boundaries 

2/1 & y2, 

• where q is the charge of the quark on such a line. Values < Sq^ >=< {q— < q >)^ = 0.22 • • • 0.25 
are obtained, and 

• where ncommon lines counts the number of quark lines intersecting both borders, as illustrated in the 
simple example given in figure 5: 

Most observables of charge fluctuations can be expressed using this basic correlation. Our fluctua- 
tion of the charges within a [— ^max. , +2/max.] box contains a combination of three such correlations. 
A simple summation yields: 

< SQ[hox]^ >= niincs entering box < < 9 >)^ (9) 

where nunes entering box is the number of quark lines entering the box. 

7 Charge Fluctuations in Equilibrium Models 

Let us use this relation to consider the prediction in more detail. In the thermalized limit with an 
infinite reservoir outside and a finite number of quarks inside, all quark inside will connect with 
quark lines to the outside as shown in figure 6. The dispersion of the charge transfer is therefore 
proportional to the total number of quarks or particles inside. 

In an "hadron gas" all particles contain two independent quarks each contributing to the 
fluctuation with roughly 1/4 yielding the estimate 1/2 as required by equation 2 . In the factorizing 
limit mesons have a 50% chance to be charged; possible baryon contribution require only a minor 
correction. 

For the " quark gluon gas" ignoring hadronization one obtains one quark charge fluctuation 1 /4 
for each charged parton. Equation 2 is drastically changed by a factor of 4. It is, however, not 
easy for this prediction to survive hadronization. If hadronization would just group initial partons 
into hadrons, the factorizing hadron gas description would stay completely unchanged. For the 
reduction it is essential to have only a single quark line contributing to the fluctuation. Only one 



5 





-11 <L'' q -i - 


Pi « 










: — * 












- ' 9 



max max 



Figure 6: Quark lines entering the box in the thermodynamic limit 
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Figure 7: Quark lines entering the box with local compensation of charge 



quark of each hadron has to originate in the primary partonic process and the other quark has to 
originate in local fluctuations and has not to contribute. The mechanism requires a sufficiently 
large box so that short range correlation can be avoided. 



8 The expanding Box 

Let us first consider the limit of a tiny box. Looking only at the first order in Ay one trivially 
obtains in any model 

<SQ^> / < TVchargcd >= 1 (10) 

which corresponds to the hadron gas value. If the box size increases to one or two units of rapidity 
on each side this ratio will typically decrease, as most models contain a short range component 
in the charge fluctuations. The decreasing is not very distinctive. In hadron hadron scattering 
processes such short range correlations are known to play a significant role and there is no reason 
not to expect such correlations for the heavy ion case. 

After a box size passed the short range the decisive region starts. In all global equilibrium 
models jl^, |o[ ^ the ratio will have to reach now a flat value. The only correction comes from 
overall charge conservation. If the box involves a signiflcant part of the total rapidity, it will force 
the ratio to drop by a correction factor 





factor = / PchZgedy / / PchZgcdy Ocl - ymax./^kin.max. (H) 



9 Charge Fluctuations in String Models 

This flatness is not expected in string models and numerical calculations indicate a manifestly 
different behavior. Only quark lines which intersect boundaries and which contribute to the 
charge measure have to be considered. String models contain local compensation of charge. Only 
contributions of lines originating around the boundaries (as illustrated in the figure 7) will appear. 
If the distance is larger than the range of charge compensation the dispersion will no longer 
increase with the box size. The total contribution will now be just proportional to the density of 
the particles at the boundaries: 

< 6Q^ > CX /5chargcd(ymax.)- (12) 

It now just counts the number of strings. 
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Figure 8: Comparison of the dispersion of the charge distribution with the density on the boundary 
of the considered box for central gold gold resp. lead lead scattering at RHIC and LHC energies. 



This resulting scaHng is illustrated in a comparison between both quantities in (12). Shown 
in figure 8 are the predictions of the Dual Parton model implementation DPMJET ||2^ for RHIC 
and LHC energies. The agreement is comparable to the proton-proton case shown in figure 9. The 
proportionality is expected to hold for a gap with roughly ^6y > 1 as for smaller boxes some of the 
quark fines intersect both boundaries. For large rapidity sizes there is a minor increase from the 
leading charge flow Ql originating in the incoming particles. In a more careful consideration ||^ 
one can subtract this contribution 

leading charge migration — <Ql> (1- < Ql >) (13) 

and concentrate truly on the fluctuation. 

A rough estimate of the relative size - with a width of neighboring string break ups and a width 
from resonance decays - leads to consistent values |2^ . 

10 Bleicher, Jeon, Koch's observation 

In a recent publication Bleicher, Jeon, Koch showed: 

• The overall charge conservation cannot be ignored at SPS energies 

• It obfiviates in this energy range the distinction between even the most extreme models 
including string models and statistical models with hadronic equilibrium. 

They showed that their string model prediction]^ coincides with the expectation of a statistical 
model of hadrons. Our string model DPMJET supports this conclusion for the SPS energy range 
as it also obtains fluctuations consistend with the "statistical" expectation. 

^ In the energy above ^/s = 5 GeV their UrQMD code is described Isil to be dominated by string fragmentation. 
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Figure 9: Comparison of the dispersion of the charge distribution with the density on the boundary 
of the considered box for proton-proton scattering at RHIC and LHC energies. 

While forward-backward hemisphere charge fluctuations were meaningful in the FNAL-SPS 
energy region, the fluctuations of charges into a central box contain two borders and require a 
correspondingly doubled rapidity range. Unfortunately this means a lot in energies. They are not 
available at SPS energies. 

11 A reference model with statistical fluctuation 

It was argued that the experimental results should be "purified" to account for charge con- 
servation. We basically aggree with such a correction. Given the uncertainites the correction 
obviously has to stay on the modeUing side. 

To obtain such a reference model we randomize charges a posteriori. To accurately conserve 
energy and momentum it was done separately for pions, kaons and nucleons^. Using DTMJET 
events for RHIC and LHC energies for proton-proton and central lead-lead collisions we obtain 
the "statistical" prediction shown in figure 10. 

To check consistency we employed the correction factor proposed by 



and indeed obtained the fiat distribution with the expected "hadron gas" value. 
12 String model versus randomized "hadron gas" 

Taking the DPMJET string model and the randomized "hadron gas" version as extreme cases (with 
the "parton gas" somewhere in between) we can investigate the decisive power of the measure. As 
shown in figure 11 we find that there is a measurable distinction at RHIC energies and sizable at 
LHC energies. 

Obviously, the method can also be directly applied to experimental data, at least in a simplified way. 
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Figure 10: Charge fluctuations with a posteriori randomized charges for p-p scattering and the 
most central 5% in Pb-Pb scattering at RHIC energies {^/s = 200 A GeV) and at LHC energies 
{y^ = 6000 A GeV) . The results are also shown with a correction factor to account for the overall 
charge conservation. 




Figure 11: Comparison of the charge fluctuations obtained in a string model DPMJET with a 
model using a posteriori randomized charges for p-p scattering and the most central 5% in Pb-Pb 
scattering at RHIC energies {y/s = 200 A GeV) and at LHC energies {y/s = 6000 A GeV). 
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Figure 12: The b dependence of the charge fluctuations obtained in a string model DPMJET for 
p-p scattering and the most central 5% in Pb-Pb scattering at RHIC energies {^/s = 200 A GeV). 



The similarity of p-p and Pb-Pb scattering is not surprising. The distinction between both 
cases is expected from the difference in collective effects. The data for p-p scattering are known to 
follow the string models, while interaction of comovers, or medium range or complete equilibrium 
will move the curve upward to a more statistical situation. These effects are presently outside of the 
model. A measured charge correlation between both extremes will directly reflect the underlying 
new physics. 

13 The b dependence of the charge fluctuations 

A similar result is obtained when the dependence on the centrality is studied. Without collective 
effects no such dependence is expected as observed in the model calculation shown in flgure 12 {b is 
the impact parameter). This experimentally measurable centrality dependence allows to directly 
observe collective effects without reference to model calculations and underlying concepts. 

Conclusion 

In the paper we demonstrated that the dispersion of the charge distribution in a central box of 
varying size is an extremely powerful measure. 

Within the string model calculation the dispersion seen in relation to the spectra shows no 
difference between simple proton-proton scattering and central lead lead scattering even though 
both quantities change roughly by a factor of 400. 

The dispersion allows to clearly distinguish between conventional string models and hadronic 
thermal models for a rapidity range which could be available at RHIC energies. In many models 
the truth is expected to lie somewhere in between and it is a reasonable expextation that the 
situation can be positioned in a quantitative way. 



10 



References 

U. Idschok et al. [Bonn-Hamburg-Munich Collaboration], Nucl. Phys. B67, 93 (1973). 
A. Bialas, K. Fialkowski, M. Jezabek and M. Zielinski, Acta Phys. Polon. B 6, 39 (1975). 
J. Whitmore, Phys. Rept. 27, 187 (1976). 
L. Foa, Phys. Rept. 22, 1 (1975). 

R. Brandelik et al. [TASSO Collaboration], Phys. Lett. B 100, 357 (1981). 
C. Quigg and G. H. Thomas, Phys. Rev. D7,2757 (1973). 
C. Quigg, ' Phys. Rev. D12, 834 (1975). 
F. W. Bopp, Riv. Nuovo Cim. 1, 1 (1978). 
S. P. Misra and B. K. Parida, Pramana20, 375 (1983). 
J. Ranft, Fortsch. Phys. 23, 467 (1975). 
T. T. Chou and C. N. Yang, Phys. Rev. D7, 1425 (1973). 
J. L. Newmeyer and D. Sivers, Phys. Rev. D 10, 204 (1974). 
K. F. Loe, K. K. Phua and S. C. Chan, Lett. Nuovo Cim. 18, 137 (1977). 
E. N. Argyres and C. S. Lam, Phys. Rev. D 16, 114 (1977). 
C. B. Chiu and K. Wang, Phys. Rev. D 13, 3045 (1976). 
J. Bias de Deus and S. Jadach, Phys. Lett. B 66, 81 (1977). 
M. Jezabek, Phys. Lett. B 67, 292 (1977). 

M. Asakawa, U. Heinz and B. Muller, Phys. Rev. Lett. 85, 2072 (2000). 



S. Jeon and V. Koch, Phys. Rev. Lett. 85, 2076 (2000) |hep-ph/000316q] . 



S. .Jeon and V. Koch, Phys. Rev. Lett. 83, 5435 (1999) ||nucl-th/9906074|] . 



M. Bleicher, S. Jeon and V. Koch, Phys. Re v. C62, 061902 (2000) |[hep-ph/0006201j] ; V. Koch, 
M. Bleicher and S. Jeon, |nucl-th/0103084 

F. W. Bopp,Nucl. Phys. B191, 75 (1981). 



A. Bialas, arXiv:hep-ph/0203047. 



K. Fialkowski and R. Wit,liep-ph/0006023, "Charge fluctuations in a flnal state with 
QGP," [hep-ph/0101258| . 

J. Ranft Phys. Rev. D 5 10 64 (1995); J. Ranft, |hep-ph/9911213| (Siegen preprint SI-99-5); 
J. Ranft, hep-ph/9911232| (Siegen preprint SI-99-6). 

Figure 2 of initial version: F. W. Bopp and J. Ranft, |hep-ph/0009327 . 

K. Adcox [PHENIX Collaboration], |arXiv:nucl-ex/0203014 . 



F. W. Bopp and J. Ranft, Eur. Phys. J. C 22, 171 (2001) |arXiv:hep-ph/0105192 
R. Baier and F. W. Bopp, Nucl. Phys. B79, 344(1974). 
P. Aurenche and F. W. Bopp, Nucl. Phys. B119, 157 (1977). 
M. Bleicher et al., J. Phys. G25 (1999) 1859 [ ►iep-ph/9909407 |. 



11 



